HESSENBERG DECOMPOSITION OF MATRIX FIELDS 
AND BOUNDED OPERATOR FIELDS 
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Abstract. Hessenberg decomposition is the basic tool used in computational linear algebra 
to approximate the eigenvalues of a matrix. In this article, we generalize Hessenberg decompo- 
sition to continuous matrix fields over topological spaces. This works in great generality: the 
space is only required to be normal and to have finite covering dimension. As applications, we 
derive some new structure results on self-adjoint matrix fields, we establish some eigenvalue 
separation results, and we generalize to all finite-dimensional normal spaces a classical result on 
trivial summands of vector bundles. Finally, we develop a variant of Hessenberg decomposition 
for fields of bounded operators on a separable, infinite-dimensional Hilbert space. 
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Introduction 

A Hessenberg matrix is a complex square matrix that is zero outside of the upper triangle 
and first subdiagonal: in other words, it is a matrix of the form 




Some authors [2j call a Hessenberg matrix unreduced if all the coefficients on the first subdiag- 
onal are nonzero. In this article, most of our Hessenberg matrices will happen to have positive 
(nonzero) coefficients on the first subdiagonal. 

A Hessenberg decomposition of a square matrix m is a decomposition 

m = uhu* 
1 



2 



BENOIT JACOB 



where u is unitary and h is Hessenberg. This decomposition is very important to computa- 
tional linear algebra, because it is the main step toward unitary triangularization, which is 
decomposing a matrix m as 

m = vtv* 

where v is unitary and t is triangular. Indeed, for nxn matrices, the Hessenberg decomposition 
can be computed exactly in O(n^) operations, and from there, to any given degree of precision, 
the unitary triangularization can be computed approximately in just O(n^) operations, and 
even just 0{n) operations in the self-adjoint case. This is explained in [2j, Chapters 7 and 8. 
Unitary triangularization is of course very important as it gives the eigenvalues, and amounts 
in the self-adjoint case to diagonalization. 

The basic observation of this article is that Hessenberg decomposition can be generalized to 
work for matrix fields, that is, continuously over a topological space, and that that can be used 
to derive new structure results on matrix fields. Our main Hessenberg decomposition results 
are summarized in Theorem 13.11 

Theorem (See Theorem 13. ip . Let ^ d < oo. Define c as: 

f ifdi^l 
c = < 2 i/2^(i^3 
( ifd^A. 

Let X be a normal space of covering dimension d. Let n ^ 1. Let f £ C(X, M„). Let 
e £ C(X, M>o). There exist g £ C(X, M„) and u £ C(X, U^) such that for all x £ X, 

• g{x) — f{x) is self-adjoint, 

• \\f{x)-9{x)\\ < eix), 

• the matrix {ugu*){x) belongs to H^~'^ (see Definition \2. 
We then derive some applications of that theorem. 

First, Theorem 14.11 is a general decomposition result for self-adjoint matrix fields: it says 
that if we allow an arbitrarily small perturbation, conjugation by a unitary field and a rank 
one perturbation, then any self-adjoint matrix field decomposes as 

Aipi H h XkPk + r 

where the pi are mutually orthogonal rank one projections, the Aj are real-valued functions, 
and r is orthogonal to the pi and is zero outside of a block of size roughly half the covering 
dimension of the base space. 

Corollary 14.31 gives a different kind of structure result in the special case of base spaces of 
covering dimension 3. 

Then, we establish some eigenvalue separation results: a generic result as Theorem 15.11 fol- 
lowed by specialized theorems in low dimensions: Theorem 15.31 shows that over a base space 
of dimension at most 2, the eigenvalues can be completely separated, generalizing a result of 
Choi and Elliott [1, Theorem 1]; and Theorem [53] gives an optimal eigenvalue separation result 
when the base space has dimension 3 or 4 — indeed, as we will recall in section \5\ complete 
separation of eigenvalues is not possible in general as soon as the dimension of the base space 
is more than 2. 

We then establish a specialized variant of Hessenberg decomposition for the case of projec- 
tion fields: that is Theorem 16. 2[ As an application (Theorem 17. ip . we generalize a classical 
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theorem on vector bundles according to which over a space of dimension d, any complex vector 
bundle of rank n has a trivial summand of rank roughly n — d/2. This has long been known 
for CW-complexes, and the case of compact Hausdorff spaces can be reduced to that case, as 
noted by Phillips [6, Proposition 4.2]. Our Theorem 17.11 shows that this actually works for all 
normal spaces. 

In the last section of this article, we replace matrices by bounded operators on an infinite- 
dimensional, separable Hilbert space. Here, a classical analogue of the notion of being diag- 
onalizable with eigenvalues of multiplicity 1, is the notion of a cyclic operator. An operator 
is said to be cyclic if there exists a cyclic vector for it. It is well-known (see [3], Chapter 18) 
that cyclic operators do not form a norm-dense subset of bounded operators, and that their 
complement is norm-dense. Therefore, contrary to the case of finite matrices, here, already in 
the case where X is a point, it is nontrivial to say anything. Fortunately, there still is a notion 
of Hessenberg operators (see Definition 18. ip . that have many interesting properties. They are 
cyclic, and it is well-known (Lemma 18. 4p that cyclic operators are exactly the operators of the 
form 

u*hu 

with u unitary and h Hessenberg. Moreover, the self-adjoint Hessenberg operators are exactly 
the Jacobi operators, whose spectral theory and inverse spectral theory are so well-developed 
(see [7])- We obtain the following result: 

Theorem (See Theorem lS.Sp . Let X be either a compact space, or a finite- dimensional normal 
space. Let H = ^^(N*). Let f be a strongly continuous map from X to B{H). For any 
£ G C(X, M>o), there exist maps g and v from X to B[H) such that, letting h = v*gv, the 
following properties are satisfied: 

• v{x) is an isometry for all x £ X . 

• h{x) is Hessenberg for all x £ X . 

• V, g, h are strongly continuous. 

• ||/(x) — ^(x)!! < e(x) for all x £ X. Moreover, f — g is norm- continuous, compact, and 
self-adjoint. 

Notice that in this theorem, v{x) is only an isometry, not a unitary in general. The statement 
would be wrong if we said "unitary" instead of "isometry" , again because cyclic operators do 
not form a norm-dense subset of bounded operators. 

The author thanks George Elliott and Leonel Robert for helpful conversations, and Etienne 
Blanchard and N. Christopher Phillips for helpful comments on an early version of the present 
article. 



Notations and terminology 

N is the set of non-negative integers, N* is the set of positive integers, IR>o is the set of 
positive real numbers. For x G M we let [x] denote the smallest integer n such that n ^ x. 

C(X, Y) is the set of all continuous maps from a space X to a space Y. The notion of 
topological dimension used throughout this article is the Lebesgue covering dimension, whose 
definition is recalled in Section [H By the dimension of a space X, we mean its covering 
dimension, which we denote by dimX (see section [T]). 

M„ is the set of complex n x n matrices. The notation mij means the (i,j)-th coefficient of 
the matrix m, that is the coefficient at row i and column j, with the numbering starting at 1. 
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M^'^ is the subset of self-adjoint matrices, U„, is the subset of unitary matrices. We will 
introduce a notation in Definition 12.11 and a notation BH^ in Definition 16. li 

A matrix is called positive if it is of the form aa* for some matrix a. This is what is 
sometimes called "positive indefinite". Thus we consider a positive matrix, even though we 
do not consider it a positive number, an inconsistent by usual terminology. A matrix m is 
negative if —m is positive. 

Given two matrices a G M„ and b G Mp, we let a © 6 denote the block-diagonal matrix 
(q ^) G M„_(.p. We let M„ © Mp denote the subalgebra of M„+p consisting of all matrices of 
that form. By convention, we let Mq = and M„ © Mq = M„. 

On a Hilbert space H, we only use the norm topology. The scalar product is denoted (•, •) 
and is linear in the first variable. We let B{H) denote the set of bounded operators on H. On 
B(H), by the strong and weak topologies we mean the operator topologies, i.e. the SOT and 
WOT. The adverbs strongly and weakly refer to these topologies. For example, to say that 
a map / : X — > B{H) is strongly continuous means that for all ^ £ H, the map x i— > /(x)^ 
is continuous. When we apply operator terminology to maps f : X ^ B(H), we mean it 
pointwise. For example, to say that / is compact means that f{x) is a compact operator for 
all X € X. 

In section [8] we will fix ourselves H = ^^(N*), we will let (ej)jgN* denote its standard Hilbert 
basis, and we will let Oij = {aci, ej) denote the ij-th. matrix coefficient of an operator a G B{H). 

1. Lemmas in topological dimension theory 

Let us first recall the classical notion of Lebesgue covering dimension, often called topological 
dimension. Given a topological space X, a refinement of an open covering {Ui)i^i of X is an 
open covering of X such that for all j £ J, there exists i £ I such that Vj C C/j. 

Definition 1.1 (See [4J). Let d G N. A topological space X is said to have dimension at 
most d if any open covering of it has a refinement {Vj)j^j such that for all x £ X , the set 
{j £ J, X £ Vj} has at most d + 1 elements. 

Obviously, dimX is then defined as the smallest d such that X has dimension at most d, or 
oo if no such d exists. It is true that dimM*^ = d. 

Let us also recall the notion of a normal space, which is the only "separation axiom" that 
we will use throughout this article: 

Definition 1.2. A topological space is said to be normal if any two disjoint closed subsets have 
disjoint neighborhoods. In other words: for all closed subsets F,G of X, if F Ci G = then 
there exist open subsets U,V of X such that F C U , G <ZV , U nV = . 

Recall the following classical theorem in dimension theory: 

Theorem 1.3 (See [1], Theorem VII. 9). Let X be a normal space. Let n G N. The following 
are equivalent: 

(1) dimX ^ n. 

(2) For any f £ C{X, [0;!]"+^), for any e > 0, for any y £ [0; there exists g £ 
C{X, [0; 1]"+^) such that \\f - g\\ < e and y g{X). 

We will need some variants and refinements of the implication [T]=>[2] in the above theorem. 
In order to obtain them, we will simply adapt the classical proof of that theorem. The main 
technical lemmas used in that proof are the following: 

Lemma 1.4 (Urysohn's Lemma). Let X be a normal space. Let F, G be disjoint closed subsets 
of X . There exists a continuous function ip £ C{X, [0; 1]) such that <f{F) = and (p{G) = 1. 
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Here, by i^iF) = A we mean that (p{x) = A for all x £ F. 

Lemma 1.5 (See [4J, VII. 4. B). Let n G N. Let X be a normal space such that dimX ^ n. 
Let Ui, . . . , Un+i be open subsets of X. Let Fi, . . . , Fn+i be closed subsets of X. Suppose that 
Fi C Ui for all i. It follows that there exist open subsets Vi, . . . , Vn+i and Wi, . . . , Wn+i of X 
such that 

FiCViCViCWiCU fori = I,..., n + 1 

and 

n+l 

fl {Wi - Vi) = 0. 

1=1 

Having recalled these classical lemmas, we can now start proving the lemmas that we will 
need. The proof of the following lemma follows very closely the classical proof of Theorem 1 1.3[ 

Lemma 1.6. Let n G N. Let X be a normal space such that dimX ^ n. For any f £ 

C(X,M"+i) and for any e G C(X,M>o), there exists g G C(X,M"+i) such that for all x £ X, 
\\g{x) — f{x)\\ < e{x) and g{x) ^ 0. 

Proof. Write / = (/i, . . . , fn+i) where the fi G C(X, M) are continuous functions. For 1 ^ i ^ 
n+l, let 

Fi = {x G X; f,{x) ^ e(x)} 
Gi = {x G X; f,{x) ^ -e{x)} 



Since Fi C X — Gi for all i, it follows from Lemma [1.5 1 that there exist open subsets Vi, . . . , Vn+i 
and Wi, . . . , Wn+i of X such that 

FiCViCViCWiC X -Gi for i = 1, . . . , n + 1 

and 

n+l 

(1) n - ■^o = 0- 

i=l 

Since X is normal, by Urysohn's Lemma 11.41 for all i there exists a continuous function ipi : 
X [—1; 1] such that ipi{Vi) = 1 and ifi{X — Wi) = —1. Since Fi C Vi, we have fi{Fi) = 1. 
Since Gi C X — Wi, we have (pi{Gi) = —1. We may therefore define a continuous function 
gi : X ^ W hy letting, for all x £ X, 

( \ ^ i fii^) if X G FiUGi 

^'^^^ " \ e{x)^i{x) iix^FiUGi 

We now define g : X ^ M"+^ by letting g{x) = {gi{x), . . . ,gn+i{x)). It is clear that ||(?i(x) — 
/j(x)|| ^ 2e(x) for all x and all i, and therefore 

\\g{x) - /(x)|| ^ 2Vn + le{x) for ah x G X 

It remains to show that g{x) ^ for all x £ X. Suppose that g{x) = for some x £ X. Then 
gi{x) = for i = 1, . . . , n + 1. It follows that x U Gi, so that gi{x) = e{x)(pi{x) = 0. 
Since e(x) > 0, it follows that aj(x) = 0. This in turn entails that x ^Vi and x ^ X — Wi. 
Therefore, x £ Wi — Vi for alH = 1, . . . , n + 1, contradicting equation ([T|). □ 

Here is now a variant where instead of avoiding just one point, we now avoid any finite 
number of maps pointwise. This is the first, but not the only place where it is useful to have 
introduced the non-constant e in Lemma 11.61 This variant will be especially useful when we 
will establish separation of eigenvalues over low-dimensional spaces (Theorems 15.31 and I5.4p . 
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Lemma 1.7. Let n G N. Let X be a normal space such that dimX ^ n. For any f £ 

C(X,R"+i), for any k E N*, for any hi,...,hk G C(X,M"+i), and for any e G C(X,M>o), 
there exists g G C(X, M""*"^) such that for all x G X, \\g{x) — /(x)|| < e{x) and 

g{x) {hi{x), . . .,hk{x)}. 

Proof. Let us work by induction on k. The case k = 1 obviously reduces to Lemma 11.61 by just 
translating by —hi. Let us now suppose the result to hold for a fixed k and let us establish it 
for k + 1. Applying the result for that k and for e/2, we obtain a map 7 G C(X, R"+^) such 
that for ah x e X, \\j{x) - f{x)\\ < e{x)/2 and 7(x) . . . , hk{x)]. Now let 



e'{x) = min ( e{x)/2, min ||7(x) — ) 

y i=l,...,k J 



Applying the result for k' = 1, for h[ = h^+i, and for e' , to the map 7, we obtain a new map 
g G C(X,M"+^) such that for all xeX, 

\\g{x)-f{x)\\ <e{x)/2 + e{xy ^e{x) 

and g{x) 7^ hk+i{x), and such that moreover \\g{x) — j{x)\\ < \\g{x) — hi{x)\\ for i = 1, . . . , fc, 
which entails that g{x) / hi{x) iov i = 1, . . . , k. □ 

Here is another variant for maps into topological manifolds, which we will have to apply to 
the spheres 5" in our main Hessenberg reduction process (Theorem 13. 2p . 

Lemma 1.8. Let n G N. Let X be a normal space such that dimX ^ n. Let Y be a metric 
topological manifold of dimension n+ 1, with metric denoted by d. Suppose that Y has an atlas 
consisting of bi-Lipschitz charts. Let Z be a discrete subset ofY. For any f G C{X,Y) and for 
any e G C(X, M>o), there exists g G C{X,Y) such that for all x £ X, d{g{x), f{x)) < e{x) and 
g{x) Z. Moreover, for any neighborhood W of Z inY , g may be chosen so that g{x) = f{x) 
for allx^ f'^iW). 

Proof. It is enough to prove that for any z £ Z, for any neighborhood U of 2; in Y, there is an 
open subset V of U such that V Ci Z = {z}, and a continuous map g G C{X,Y) such that for 
aU X e X, d{g{x), f{x)) < e{x), g{x) Z, and if x f~^{V) then g{x) = f{x). 

So lei z £ Z and let [/ be a neighborhood of z in Y . By assumption, there exists an open 
subset V oiU such that V r\ Z = {z} and a bi-Lipschitz map ip -.V ^ ioi some open subset 
of M"+^. Moreover we may choose ip and 17 so that 'p{z) = 0. For any r > 0, let denote 
the open ball in M""*"^ of radius r centered at in IR"'^^. Choose > so that 

Bzrf C Vt. 

Let A = f~^{y). Let /' G C(A, 17) be defined as f = ip o f . By Lemma [1.6| there exists a map 
h' G C(^,M"+i) such that \\h' - f'\\ < r? and h'{A). Now define a map g' G C{A, B2rj) as 
follows. Let X G ^. If f'{x) G Br,, then set g'{x) = h'{x). If f'{x) B2r,, then set g'{x) = fix). 
Otherwise, we have r] ^ ^ 2ri, we let 

V 

and set g'{x) = tf'{x) + (1 - t)h'{x). Define g G C{A,V) as g = ip~'^ o g' . Notice that g 
agrees with / outside of A, hence g can be extended to all of X by letting g(x) = f{x) for all 
A. □ 
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2. A CLASS OF HeSSENBERG-LIKE MATRICES 

Definition 2.1. For n G N* and /c G {0, . . . , n}, let be the set of all matrices m G M„ such 
that 

iTLij = whenever j ^ k and i ^ j + 2 

and 

ruij G M>o whenever j ^ k and i = j + 1. 

Notice that = M„ and that H^~^ = is the set of ah matrices in M„ that are Hessenberg 
and that have positive coefficients on the first subdiagonal. In particular, since the subdiagonal 
coefficients are nonzero, these are called unreduced Hessenberg matrices. More generally, for 
any k G {0, . . . the matrices in are "unreduced Hessenberg in their k first columns". 
They are the matrices of the form 




* * * 



where a "+" indicates a coefficient in M>o. 

Let us make the following observations on the eigenvalues of matrices in H^. 

Proposition 2.2. Let n G N* and k G {0, . . . ,n}. Let p = max(n — A;, 1). Let x G H^. It 
follows that all eigenvalues of x have multiplicity at most p. 

Proof. Let A be an eigenvalue of x. Notice that x—X 1„, still belongs to H^. Let q = min(n— 1, k). 
It follows from the definition of that the q first columns of x — A 1„ are linearly independent, 
hence rr — A 1„ has rank at least q, hence dim ker (x — Xln) ^ n — q = p, hence A is of multiplicity 
at most p. □ 

The proof of our second lemma closely follows the proof of a theorem well-known as the 
Sturm Sequence Property (see [2], Theorem 8.5.1). 

Proposition 2.3. Let n G N* and k G {0, . . . ,n}. Let p = max(n — k — 1,1). Let x G be 
self-adjoint. It follows that at most p eigenvalues of x have multiplicity more than 1. 

Proof. The statement is trivial if n ^ 2, so let us assume that n ^ 3. Also the statement 
is trivial if /c = 0, so let us assume that k ^ 1. For i G {0, ... ,n — 1}, let Cj denote the 
bottom-right {n — i + 1) x {n — i + 1) corner of x. Notice that ci = x and that c„ is the 1x1 
matrix [Xnn 

). Let 

Pi{X) = det(ci - A In-i+i) 

be its characteristic polynomial. Expanding determinants with respect to the first column, one 
easily establishes the relation 

(2) pi{X) = XiiPi+i{X) - Xi+i jPi+2(A) for 1 ^ z ^ min(A;, n - 2) 
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It follows (see [9], §47) from the minimax principle that if Ai ^ . . . ^ A„ are the eigenvalues of 
X = ci and /Ui ^ . . . ^ fJ-n-i are the eigenvalues of C2, then the following Interlacing Property 
holds: 

(3) Ai ^ /il ^ A2 ^ . . . ^ fln-l ^ A„. 

Now let A be an eigenvalue of x, and suppose that A has multiplicity more than 1. It follows 
from the inequalities ([3]) that A is not only an eigenvalue of x = ci but also an eigenvalue of 
C2- In other words: 

Pi(A) =p2(A) = 0. 

It follows from equation ([2|) that X2 iP3(A) = 0. Since x G and /c ^ 1, we have 2:2,1 7^ and 
it follows that 

P3(A) = 0. 

We can continue applying equation ([2]) iteratively as long as i ^ min(A;,n — 2). So, letting 
j = min(A:, n — 2), in the end we obtain 

= Pi(A) = P2(A) = • • • = Pj+2(A). 

Thus, A is an eigenvalue of Cj+2, which is a square matrix of size n — j — 1 and therefore can't 
have more than n — j — 1 = max(n — k — 1,1) distinct eigenvalues. □ 

3. Hessenberg reduction of matrix fields 

The goal of this section is to establish the following theorem: 

Theorem 3.1. Let ^ d < 00. Let c be defined as follows: 

ifdi^l 
2 i/ 2 ^ d ^ 3 

ifd^A. 

Let X be a normal space of covering dimension d. Let n ^ 1. Let f G C(X, M„). Let 
e G C(X, M>o). There exist g £ C(X, M„) and u G C(X, U„) such that for all x £ X, 

• g{x) — f{x) is self-adjoint, 

• \\f{x)-g{x)\\ < e{x), 

• the matrix {ugu*){x) belongs to H"~'^ (see Definition \2. 

We will establish it as three separate propositions: 13.21 13. 3|, 13. 4[ 

Proposition 3.2. Let n G N* and d G N. Let k = n — \d/2\ — 1. Let X be a normal space 
such that dimX ^ d. For any f G C(X, M„) and e G C(X,M>o), there exists g G C(X, H^) 
and u G C{X, U„) such that for all x £ X , (/ — u*gu){x) is self-adjoint and has norm less than 
e{x). 

Proof. First notice that it is sufficient to prove that for all p G {0, . . . , A; — 1}, / G C(X, H^), 
and e G C(X, M>o), there exists g G C(X, H^"''^) and u G C(X, Un) such that for all x £ X, 
(/ — u*gu){x) is self-adjoint and has norm less than e{x). 

Let b be the block inside the (p+ l)-th column of / starting at the (p + 2)-th row and ending 
at the n-th row. In other words, 

b= : 

V bn,p+l J 
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In the following diagram, we represent the matrix / G M„(C(X)) with the block b framed 
inside it: 

(4) 




* * 



In this diagram, the p first columns are represented ending with zeros. Of course, this only is 
an accurate representation if p ^ 1. In the case p = 0, one should simply imagine that the 
column with the frame is the first column. 



Notice that since every coefficient fij is a function in C{X), the block b may be seen as a 
continuous map 

Since p^k — l^n — d/2 — 2, we have 

(5) d^2n-2p-A. 

By Lemma 11.61 and inequality there exists a continuous map 
for all X e X, \\b{x) - b'{x)\\ < e{x)/n and b'{x) / 0. 

For X G X, let r{x) = \\b'{x)\\ £ IR>o, and let = b'{x)/r{x). Notice that b'{x) = r{x)(3{x) 
and that (5{x) belongs to the unit sphere 5'2n-2p-3 -^^ ]^2(n-p-i)_ 

By Lemma [L8l and inequality ([5]), there exists a continuous map (3' : X ^ ^2n-2p-3 g^^i^ 
that for all X G X, 

eix) 



?'(x)-/3(x)|| < 



nr{x) 



and 



13' {x) + 







p2n-2p-2 



Now let us see S*^" '^^ ^ as the unit sphere in C" ^ so that is a map into C" ^ ^. 
Define a map b" : X ^ (^n-p-i |-,y letting 6"(x) = r(x)/3'(x). Notice that for all x G X, 
||6"(x) — 6(x)|| < £{x)/n and, when seeing b" as a map into IR^^""^"-'^^ 



(6) 



b"{x) / 







for all X G X and A G M<o. 
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Now let us construct a self-adjoint matrix field 5 such that adding 6 to f has the effect of 
replacing the block b by b". This is trivially done by letting, for all x £ X, 



5ij{x) = < 



b'-_„_iix) - bi-p-i{x) if i ^ p + 2 and j = p+l 



bj_p_i{x) - bj-p-i{x) if j ^ p + 2 and i = p+l 







otherwise 



Notice that \\5{x)\\ < e(x). Let 



It remains to construct a unitary field u G C{X, U„) such that {uf"u*){x) is always a matrix in 
Hn^^. We will construct u as a Householder reflection, a technique widely used in computational 
linear algebra. Let us look back at the diagram ^ applied to the matrix /". The block framed 
there is exactly b" . This is a map from X to C"^^^^, and we know by ([6]) that it avoids the 
closed real half-line generated by (—1,0, ... ,0). Therefore, we may define a continuous map 
h:X^ by letting 



h{x) 



b"{ 



\b"(x] 



+ 



/1\ 



for all X E X, 



and we then have h{x) ^ for all x £ X. The vector h{x) is what is known as a Householder 
vector. Let uq{x) G U„_p_i be the orthogonal reflection around the complex line generated by 
h{x). We have 

■ \h"{x) 




UQ{x)b"{x) 



for all X G X. 



V / 

Moreover, the map : X ^ U„_p_i thus defined is continuous. Finally we obtain the wanted 
unitary field u G C(X, U„) by letting, for all x G X, 



n(x) 



P+l 







uq{x)) ' 

The wanted matrix field g is then obtained as g = uf"u* 



□ 



Now that we have obtained our general Hessenberg decomposition result for spaces of ar- 
bitrary finite dimension, let us sharpen it in some special low-dimensional cases. Here is first 
a refinement for spaces of dimension at most 3. The main new tool that becomes available 
in dimensions 3 and less, is the possibility to use continuous Givens rotations while in the 
above general case we could only use Householder reflections. In computational linear algebra. 
Householder reflections and Givens rotations are the two basic types of unitaries used in almost 
all unitary reduction processes. For further details on these matters, see [2]. 

Proposition 3.3. Let n G N*. Let X be a normal space such that dimX ^ 3. For any 
f G C(X, M„) and e G C(X,]R>o), there exists g G C(X, H^"^) and u G C(X, U„) such that for 
all X £ X, (f — u*gu){x) is self-adjoint and has norm less than e(x). 
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Proof. By Theorem 13.21 we can get all what we need except that g{x) is in HJ^ and we want 
it to be in H^J-^. Since g e C(X,H;^-3), it has the form 



* 
+ 




0- 



* * * * 

+ * * * 








where "+" again denotes a coefficient that is in IR>o for aW x £ X. In the above diagram we 
framed the 2x1 block consisting of gn-i,n~2 and gn,n~2- Call it h. Thus 6 is a continuous map 
from X to C^. Since dimX ^ 3 < dimC^, it follows from Lemma 11.61 that an arbitrarily small 
perturbation on the matrix field g ensures that h{x) ^ for all x £ X. This perturbation can 
be made to be self-adjoint as in the proof of Theorem 13. 2t just apply the adjoint perturbation. 
So we assume without loss of generality that h{x) ^ for all x G X. Let r{x) = ||6(x)||. Define 
a 2 X 2 unitary field uq by letting 



uq{x) 



1 / hi{x) b2{x) 
-hix) bi{x) 



r{x) 



Notice that for all x S X we have 



uo{x)b{x) 



r{x] 
' 



ln-2 

uo{x)J' 



It follows that if we define a n x n unitary field u, usually called a Givens rotation, by letting 

u{x) 

then the matrix field ugu* has the wanted form 

+ 




* * * * 
+ * * * 








* * 



□ 



Our next Hessenberg decomposition result is for the special case of spaces of dimension at 
most 1. Here, we can easily obtain the actual unreduced Hessenberg form, i.e. our matrices are 
Hessenberg at every point and the subdiagonal coefficients do not vanish. 

Proposition 3.4. Let n G N*. Let X be a normal space such that dimX ^ 1. For any 

f £ C(X, M„) and e G C(X,M>o), there exists g G C(X,H;:^) and u G C(X,U„) such that for 
all X £ X , [f — u*gu){x) is self-adjoint and has norm less than e{x). 
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Proof. By Theorem 13. 3^ we can get all what we need except that g{x) is in HJ^ ^ and we want 
it to be in H". Since g e C(X, H™"^)) it has the form 



+ 








* * * * 
+ * * * 
+ * * 







where "+" again denotes a coefficient that is in M>o for all x € X. In the above diagram we 
framed the 1x1 block consisting of gn,n-i- All what we need is to make it positive (nonzero). 
Since dimX ^ 1 < dimC, it follows from Lemma 11.61 that an arbitrarily small perturbation of 
the matrix field g will ensm'e that gn^n-i{x) 7^ for all x ^ X. As in the proof of Theorem 13. 2^ 
this perturbation can be chosen to be self-adjoint just by applying the adjoint perturbation. 
So we may assume without loss of generality that gn,n-i{x) / for all x G X. Now let 



C(x) 



\gn,n-l{x)\ 



gn,n-l{x) 



and define a nx n unitary field u by letting u{x) 
ugu* has the wanted form. 



diag(l, . . . ,1, Cix))- It is then obvious that 

□ 



4. Application: general structure results for self- ad joint matrix fields 

Theorem 4.1. Let d G N. Let X be a normal space of dimension d. Let c = \d/2\ + 1. 
Let n G N. Assume that n > c. Let k = n — c. Let 11^ be the set of rank one (orthogonal) 



. It follows that exists u G C{X, U„), 
Xk G C(X, M) such that the following 



projections in M„. Let f G C(X,M^'"). Let e G C(X,R>o) 
g,r...,pk G C(X,M^^),pi,...,pfc G C(X,ni) andXi,... 
conditions are satisfied: 

• \\{ufu*){x) - g{x)\\ < e{x) for all x G X ; 

• g = Xipi H h XkPk + r + q; 

• The projections pi are mutually orthogonal, and are orthogonal to r; 

• The map r is zero outside of the bottom-right cx c corner: in other words, there exists 
a map rc G C{X, M^^) such that r = 0„_c © 

• q may be chosen: 

— Either so that for all x G X, q{x) has rank exactly 1 and is positive, and has norm 
at most ^/2{\\f{x)\\ +£{x)); 

— Or so that for all x £ X , q{x) has rank exactly 1, is negative, and has norm at 
most \/2(||/(x)|| + £{x)), and so that if f is positive then r,pi, . . . ,pk are positive 
too; 

— Or so that for all x G X, q(x) has rank exactly 2, trace 0, and norm at most 
||/(x)|| + e{x), and so that if f is positive then r,pi, . . . ,pk are positive too; 

• In either case, the spectral projections of q are trivial; 

• Pi + • • • + Pk is trivial too (actually it is constant equal to Ik © Oc); 

• Xi{x) > . . . > Afc(x) for all X £ X. 

Remark 4.2. Notice that while r,pi,. . . ,pk are mutually orthogonal, q is not necessarily or- 
thogonal to any of them. Also notice that when d ^ 2, there is the stronger result that the 
eigenvalues can be completely separated, see Theorem \5.3l which provides a splitting into an 
orthogonal sum of multiples of rank one projections. Finally, notice that when d = 3, Theorem 
\3.3\ allows to obtain c = 2 instead of c = 3, which is a slight improvement. 
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Proof. Apply Theorem [32] to obtain u G C(X, U„) and g G C{X,M^j^) such that for all x e X, 
\\{ufu*){x) — g{x)\\ < e{x) and g{x) G HJJ"'^- It remains to split g as advertised. Since g{x) 
belongs to H^"'^ and is self-adjoint, it has the following form: 




The two framed coefficients gk+i,k{x) and gk^k+i{x) are of course equal: call them //(x). In the 
case where we want q of constant rank 2, let 



Ofc-i e 



n 

IJL 



c-1- 



In the case where we want q of constant rank 1 and positive, let 



In the last case where we want q of constant rank 1 and negative, let 



Jc-l- 



Either way, g — q has this form: 




Let now h be the top-left kxk block and r be the bottom-right c x c block in g — as illustrated 
on this diagram: 

■ ■ 
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Implicitly we extend r by zero to make it a n x n matrix. Notice that h{x) G and therefore, 
by Proposition [221 h{x) has k distinct eigenvalues for all x £ X. Call them Ai(a;) > . . . > Afc(x) 
and let Pi{x) be the spectral projection of h{x) corresponding to Aj(x), again implicitly extended 
by zero to make it a nxn matrix. This completes the construction, and the remaining conditions 
are then easily shown to be satisfied. □ 

Following a diff'erent approach, one might ask if other structure results become available 
when the base space X has low enough dimension. For spaces of dimension up to 2, we will 
see below that eigenvalues can be completely separated: see Theorem 15.31 This shows that the 
self-adjoint field then splits as an orthogonal sum of multiples of rank one projections, which is 
about the strongest structure result that one might hope for. However, already in dimension 3 
this isn't possible anymore, as we will discuss in Section [H Hence, the following easy corollary 
of Theorem 13.31 is perhaps worth mentioning: 

Corollary 4.3. Let X be a normal space of dimension at most 3. Let f E C(X, M^'^). Let 
e G C(X,M>o). There exists g G C(X, M^*^) and u G C(X, U„) such that we have \\{ufu*){x) - 
g{x)\\ < e{x) for all x G X, and there exists an open subset U of X such that: 

• For all X £ U , g{x) has n distinct eigenvalues 

• For all X £ X — U , g{x) has the block- diagonal form 

V(x) 
A(x) 



gix) 



where g'{x) is a {n — 1) x (n — 1) has n — 1 distinct eigenvalues, and A(x) G M. 

Proof. Apply Theorem 13.31 and construct U as the set of all x G X such that gn,n-i{x) / 0. 
The statements on distinct eigenvalues follow from Proposition 12. 2[ □ 

5. Application: separation of eigenvalues of self-adjoint matrix fields 

By separation of eigenvalues, we mean results showing that, given a matrix field, by apply- 
ing a small perturbation to it, one can ensure that at every point the matrices do not have 
eigenvalues of too high multiplicity, and do not have too many eigenvalues of multiplicity more 
than 1. 



We do not claim that the approach presented in this section is the best possible, but it is 
a quick application of our Hessenberg decomposition results, and it does give the strongest 
possible conclusion at least for self-adjoint fields over base spaces of dimension up to 4. To our 
knowledge, these results are not to be found in preexisting literature, although the following 
results should be mentioned. First, Choi and Elliott [H Theorem 1] give the special case of 
compact metrizable spaces of dimension up to 2. Second, Phillips |5l Proof of Theorem 3.3] 
gives a different but related result using another method that is probably better and which 
we believe could also be used to rederive and generalize all the results that we present in this 
section. That however would require one to first establish some more advanced variants of the 
dimension-theoretic lemmas in section [T] of the present article. 

Here is a generic eigenvalue separation result, as a straightforward corollary of our Hessenberg 
decomposition. Below we will also prove finer results for spaces of low dimension. 

Theorem 5.1. Let ^ d < oo. Let X be a normal space such that dimX = d. Let n ^ 1. Let 
f G C(X,M^'^). Let £ G C(X,R>o). There exists g G C(X,M^^) such that for all x£X, 

• \\f{x)-g{x)\\ < e{x), 

• all eigenvalues of g{x) have multiplicity at most \d/2\ + 1, 

• at most [(i/2] eigenvalues of g{x) have multiplicity more than 1. 
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Remark 5.2. Notice that in particular, this means that g[x) has at least n — d—1 eigenvalues 
of multiplicity 1. Thus, for large n and fixed d, almost all eigenvalues of g[x) have multiplicity 
1. 

Proof. Apply Theorem 13.21 to / to obtain a matrix field g G C(X, H^) with k = n — \d/2] — 1. 
Since / and g — f are self-adjoint, so is g. The wanted properties on the eigenvalues of g{x) 
are then given by Propositions 12.21 and 12.31 □ 

The next two theorems are finer variants for spaces of low dimension, giving optimal eigen- 
value separation results up to dimension 4. Let us first establish total separation of eigenvalues 
when that it possible. But when is that possible? Choi and Elliot [T, Theorem 1] show that 
that is possible for compact metrizable spaces of dimension up to 2, and provide (essentially) 
the following counterexample showing that that is not possible already on the Euclidean 3-ball 

f : ^ M2(C) 

l/l — z X + iy 
2 \x — iy 1 + z 

Notice that the restriction of / to the boundary S'^ or B"^ is the Bott projection, which is a 
nontrivial rank one projection on S"^ (the corresponding line bundle is the Hopf bundle). If / 
could be approximated by normal matrix fields with distinct eigenvalues, then an approxima- 
tion of the Bott projection would extend to B^, hence would be trivial, a contradiction. 

This shows that the following is optimal as far as total separation of eigenvalues is concerned; 
below we will also give an optimal result for dimensions 3 and 4. 

Theorem 5.3. Let X be a normal space such that dimX ^ 2. Let n ^ 1. Let f € C(X, M^*^). 
Let £ G C(X,M>o). There exist g £ C(X, M^*^) such that for all x € X, ||/(x) - g{x)\\ < e{x), 
and g{x) has n distinct eigenvalues. 

Proof. The statement is trivial if n = 1, so we assume without loss of generality that n ^ 2. 
Apply Theorem 13.31 to / to obtain a self-adjoint matrix field g G C(X, H^) with k = n — 2. 
Since g G C(X, H"~^) and is self-adjoint, it has the form 



* 


+ 0- 






-0 


+ 










0^ 


\ * 












\ + 


* 


+ 










+ 


* 


* 


0- 








* 


* 



where "+" again denotes a coefficient that is in R>o for all x £ X. The framed block, which 
we will call b, is a continuous map from X to C x M: indeed, since g is self-adjoint, its diagonal 
coefficients are real. For x £ X, let Ai(x) > ... > An-i(x) denote the eigenvalues of the 
top-left (n — 1) X (n — 1) corner of g. Notice that they are all distinct by Proposition 12.21 Since 
dim(C X M) = 3 and dimX ^ 2, we may apply Lemma 11.71 showing that an arbitrarily small 
perturbation ensures that for all a; G X, 

6(x) 0{(O,Ai(x)),...,(O,A„_i(x))}. 

Moreover, this can be achieved by a self-adjoint perturbation on g, by applying a suitable 
perturbation on the other side of the diagonal. We then have that for all x £ X, either 
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9n,n-i{x) / 0, or gn,n{x) / \{x) for all i. If gn,n-i{x) / 0, then the proof of Proposition] 
gives that all eigenvalues of g{x) have multiplicity 1. If on the other hand gn,n~i{x) = and 
gn,n{x) 7^ Aj(x), then g{x) is a block-diagonal matrix of the form a® (3 where the block a is 
of size n — 1 and has n — 1 distinct eigenvalues Ai(x), . . . , An-i(x), and the block /? is just the 
1x1 matrix {gn,nix)), and we already know that gn,n{x) 7^ Aj(x) for all i, so g{x) has n distinct 
values. □ 

In dimensions 3 and 4, as we already said, complete separation of eigenvalues can't be hoped 
for, so the following is optimal: 

Theorem 5.4. Let X he a normal space such that dimX ^ 4. Let n ^ 1. Let f G C(X, M^*^). 
Let £ G C(X,M>o). There exist g G C(X,M^^) such that for all x £ X, \\f{x) - g{x)\\ < e{x), 
and g{x) has at least n — 1 distinct eigenvalues. 

Remark 5.5. Notice that in particular, this means that g{x) has no eigenvalue of multiplicity 
more than 2, and has at most one eigenvalue of multiplicity 2. 

Proof. The statement is trivial if n ^ 2, so we assume without loss of generality that n ^ 3. 
Apply Theorem 13.21 to / to obtain a matrix field h G C(X, H^) with k = n — 3. Since 
h G C{X,Yil-^), it has the form 

■ ■ 




where "+" again denotes a coefficient that is in M>o for all x G X. Let f^i{x) ^ ^2(2;) denote 
the eigenvalues of the 2x2 bottom-right corner. Let Ai(x) > . . . > Xn-2ix) denote eigenvalues 
of the (n — 2) X (n — 2) top-left corner; notice that the Aj(x) are distinct by Proposition 12.21 as 
that corner belongs to H^l2- As illustrated in the above diagram, let a{x) = hn-i,n-2{x) and 



h{x) = hn^n-2{x)- Since dimC 
that there exists a map c = 



(7) 



c{x) 




5 and dimX ^ 4, we may apply Lemma 11.71 showing 



C X C X M such that for all x G X, 



C X 



< e(x) and 



— a(x) 
-b{x) 
(Ai - ^i)(x)^ 

— a(x) 
-6(x) 
(Ai - ii2){x) 



— a(x) 
-b{x) 
^(A„_2 - /ii)(x) 

— a(x) 
-6(x) 

(A„_2 - ^^2){x), 



Now use c to define a perturbation g of h as follows. For any i,j let Eij denote the matrix 
whose (z,j)-th entry is 1 and whose other entries are all 0. For x G X, let 



g{x) = h{x) 



+ Ci{x)En-^l^n-2 + 
+ C2{x)En,n-2 + 
+ Cz{x){En-l,n-l+ En,n)- 



Ci{x) En-2,n~l 
C2{x)En-2,n 



Notice that g is still continuous and self-adjoint. Now let us fix ourselves some x G X for the 
remainder of this proof. We want to prove that g{x) has at least n — 1 distinct eigenvalues. 
For brevity, let 

a' = gn~i,n-2{x) = a{x) + ci(x) 
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and 

b' = gn,n-2{x) = h{x) + C2(x). 

We have either a' 7^ 0, or h' 7^ 0, or a' = h' = 0. The rest of the proof sphts into two cases: 
First case: a' / or 5' 7^ 0. Let r = 11 11. Let uq be the following rotation: 



no = - 



r 



—I 



Let u = ln-2 © ^0- The unitary u G is again an example of a Givens rotation. The point is 
that 



a' \ I r 



^0 \b' J - Vo, 

so that the matrix m = ug{x)u* belongs to H^~^. It is thus enough to prove that for any matrix 
m S H^~^, m has at least n — 1 distinct eigenvalues. There are two cases: either mn,n-i = 
or not. If run^n-i = then the matrix m is block-diagonal with a block in H^-i followed by a 
1x1 block, so the result follows from Proposition [2T2l If, on the other hand, mn,n—i 

^ 0, then 

conjugating by an obvious unitary ensures that m„^„„i G M>o and therefore the matrix m is 
in and the result follows again from Proposition! 



Second case: a' = b' = 0. It then follows from the definition of g that g{x) is block-diagonal 
with a (n — 2) X (n — 2) block with the n — 2 distinct eigenvalues Ai(x) > . . . > A„_2(x), followed 
by a 2 X 2 block. It then follows from ([7]) and the definition of g that the eigenvalues of the 
latter 2x2 block are not equal to Aj(x) for any i, which establishes that g{x) has at least n — 1 
distinct eigenvalues. □ 

6. The case of projections 

Projections are a special case of self-adjoint elements and so our Hessenberg decomposition 
results can readily be applied to them. There is, however, more that we can say in this special 
case, and it turns out that the following Hessenb erg-like form is more useful for projections 
than the form considered so far. 

Definition 6.1. For n G N* and k G {1,. . . ,n}, let BH^ denote the set of all matrices p in 
Mn satisfying the following conditions: 

• p is a projection (p = p* = p^)■ 
• p is a block diagonal matrix of the following form: 

p G © • • • © © M/3 

where r G N, where Ui G {1, 2} for all i, and where (3 is as follows: 



P 



n — k if ar = 1 
n — k — 1 if ar = 2. 



Notice that it is implicitly required that ai + • • • + + /? = n.. 
• Outside of the last block of size /?, all matrix coefficients of p are real nonnegative. In 
other words: pij G Mjjo whenever i,j ^ n — 0. 

Theorem 6.2. Let ^ d < 00. Define c as follows: 

ifdi^l 
c = { 2 i/2^d^3 
\f\+l ifd^A. 
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Let X he a normal space of covering dimension d. Let n ^ 1. Let p G C{X, M„) be a projection. 
There exists a projection q G C(X, M„) such that: 

• q{x) belongs to BH"~'^ for all x £ X 

• There exists u £ C(X, U^) such that q = upu* . 



Proof. Let e be any number such that 



< e < 



242 



Apply Theorem 13.11 to p and e, seeing e as a constant function on X. Let u be as given by 
the theorem. Let h = upu* . For ah x £ X, h{x) is a projection and is within distance e of g{x) 
which belongs to H"~'^. Let h'{x) be the matrix defined by h'{x)ij = h{x)ij whenever \i — j\ ^ 1 
and i,j^n — c, and h'{x)ij = otherwise. Thus h' has this shape: 




where a "+" denotes a positive coefficient. Notice that h'{x) is still within distance ne of the 
projection h{x). It follows that h'{x)'^ is within distance 4ne of h'{x). But, for all 1 ^ j ^ 
min(n — 2,n — c), the (j + 2,j)-th coefficient of h'{x)'^ is equal to /i'(x)j+ij7i'(x)j+2j+i- It 
follows that 

h' {x)j^ijh' {x)j+2,j+i < 4ne for all x G X and 1 ^ j ^ min(n — 2,n — c). 
Therefore, for all such x and j, either h'{x)j^ij or /i'(x)j+2,j+i must be smaller than 2 



ne. 

This says that for every x, the top-left (n — c) x (n — c) corner of h'{x) is "nearly" block-diagonal 
with blocks of size 1 or 2. It should be noted that this block-diagonal form (i.e. the sequence 
of block sizes) depends on x. A similar argument in the (n — c)-th column shows that for all 
X £ X, for all k ^ n — c + 2, either h'{x)n-c+i,n-c or \h' {x)i:,n-c+i\ must be smaller than 2^fn£. 
Let Q be the function defined on C by letting a(0) = and 

max (0, \ z\ — 2y/ne) 



a{z) 



for all z G C - {0}. 



Define a new matrix field h" by letting 

h"{x)ij = a{h'{x)ij) for all x G X and all i,j. 

Notice that h"{x) is within distance 2n^/2£i/2 Qf }i'(^x), hence 

\\h"{x)-h{x)\\ <6n3/V/2. 

Since e < l/{2A^n'^), it follows that \\h"{x) — h{x)\\ < 1/4. Also notice that for every x £ X, 
h"{x) has the block-diagonal shape described in Definition 16.11 for k = n — c. Again, it should 
be noted that this block-diagonal form depends on x. Anyway, for every separate x, this 
block diagonal form is preserved when we replace h"{x) by a function of it, by functional 
calculus. Since — h{x)\\ < 1/4 and h{x) is a projection, we know that 1/2 does not 

belong to the spectrum of h"{x), hence if we let X(i/2;oo) be the characteristic function of 
(1/2; oo), and let q{x) = X(i/2;oo){h" {x)) by functional calculus, it follows that q G C(X, M^) is 
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a projection, that q{x) has the same block diagonal shape as h"{x) for each x £ X, and also 
that \\q{x) - h"{x)\\ ^ \\h"{x) - h{x)\\ < 1/4 and hence 

(8) ||g(:c)-/i(x)||<i- 

The following argument is then most classical: let 

z = qh+{l- q){l - h), 

notice that it follows from inequality ([8]) that z{x) is invertible for all x G X; and construct 
the unitary v{x) of its polar decomposition as follows: 

v{x) = z{x) {z{x)* z{x))~^^'^ . 
This unitary realizes the unitary equivalence between q and h, completing the proof. □ 

7. Application: trivial summands of vector bundles 

It is a well-known principle that over a d-dimensional space, under certain conditions, com- 
plex vector bundles of rank n have trivial summands of rank roughly n — d/2. As we said in the 
introduction, this has long been known for CW-complexes, and the case of compact Hausdorff 
spaces can be reduced to that case using some advanced results of dimension theory (see [6l 
Proposition 4.2]). 

In this section, using a completely different approach, we show that this actually works for 
all normal spaces. This is obtained as a corollary of our Hessenberg decomposition results. 
Thus, this is a generalization of previously known results, with a far more elementary proof. 
By a vector bundle, we mean a locally trivial complex vector bundle. We will make use of the 
notion of finite type developed by Vaserstein ^ which, in the case of normal spaces, just means 
that there exists a finite open covering consisting of open subsets on each of which the bundle 
is trivial. 

Theorem 7.1. Let ^ d < oo. Let 7 be defined as follows: 

r ifdi^l 
j= I 1 i/2^d^3 

Let X be a normal space such that dim A" = d. Let ^ be a vector bundle of finite type over X. 
Let 

b = mindimf^. 

Suppose that 6^7 + 1. It follows that ^ has a trivial summand of rank 6 — 7. In other words, 
^ possesses cross-sections si, . . . , such that for all x G X, the family {si{x), . . . , Sb--y{x)) 
is linearly independent. 

Proof. Since ^ has finite type, by the theorems in [8], ^ is isomorphic to the column-space of a 
projection field p £ C{X, M„) for some n £ N*. By induction on b, it is enough to prove that 
if 6 ^ 7 + 1 then there exists a non- vanishing cross-section. So let us suppose that 6^7 + 1, 
and let us construct that cross-section. By Theorem 16.21 we may choose p so that 

p{x) £ BH""'^ for all x G A, with c = 7 + 1. 

Since b ^ c, for all x £ X, we have rk(p(x)) ^ c. For 1 ^ i ^ n and x £ X, let Vi{x) denote 
the i-th column of p{x), seen as a vector in C". Let (ei, . . . , e„) be the standard basis of C". 
Let us show that: 

(9) For all x G A, if ||fj(x)|| < I/V2 for all i ^ n — c, then Vn^c+i{x) = Cn^c+i- 
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Indeed, for each x G X, since p{x) belongs to BH^"*^, it decomposes as vr © p with p either 
in Mc or in Mc_i. Since p{x) is a projection, both vr and p are. Since tt is a projection, is 
block-diagonal with blocks of size at most 2, and all the columns of vr have norm strictly less 
than l/\/2, it follows that vr = 0. So p{x) = © /). Since p{x) has rank at least c, this excludes 
the case p £ Mc-i and we therefore have p € Mc. Moreover p, which is a projection of rank 
at least c, must then be the c x c identity matrix. We therefore have p{x) = On-c © Ic- This 
completes the proof of ([9]). 

For x G X, let ix denote the smallest i such that ||vj(x)|| ^ l/\/2, and let 

vix) = Vi^{x). 

Note that v is clearly a set-theoretic (i.e. perhaps not continuous) nonvanishing cross-section 
of ^. It remains to check that v really is continuous. Let xq £ X. We must show that v is 
continuous at xq. 

By Q, we know that ix^ ^ n — c + 1. Moreover, if i^o = n — c + 1 then ([9]) makes it clear 
that the mapping x ^ ix is constant on a neighborhood of xq, and hence v is continuous at xq 
in that case. Also, if \\v{x)\\ > 1/V2, the same conclusion obviously holds. 

Thus, there only remains to handle the case when ixg ^ n — c and ||u(xo)|| = l/\/2- Since 
ixo ^ n — c, the i^o'th column of p{x) falls in the area where p{x) is block diagonal with blocks 
of size 1 or 2. Since p{x) is a projection and the i^o'th column has norm l/\/2, we must be in 
presence of a diagonal block of size 2. Let j be the index such that this 2x2 block is located 
in columns j and j + 1- Thus we have 

ixo e {j,j + 1} 

and the 2x2 block in question is 

\Pix)j+i,j P(^)j+i,j+iJ ' 

This 2x2 matrix is a projection; since it has a column of norm it must be of rank 1. 

Since moreover (see Definition 16. ip its matrix coefficients are real nonnegative, it must be of 
the form 

/ a \J a — o?\ 
\\/a — 1 — a / 

with Q G [0; 1] . The norms of the two columns of such a 2 x 2 matrix are 

(10) \fa and \/\ — a 

respectively. Notice that for any value of a, at least one of the above two norms is 
This shows that there exists a neighborhood U of xq such that ^ j + 1 for all x £ U. On 
the other hand, as we already observed, the proof of ([9]) makes it clear that Vi{x) = for all 
i < j, and therefore the neighborhood U may be taken so that 

'i'x £ {j, j + 1} for all X £ U. 

Finally, since ||?;(xo)|| = l/v^, we know that one of the two numbers in (jlOp must be equal to 
l/\/2, and it follows that a = 1/2 and therefore the 2x2 diagonal block in question is equal 
to 

A/2 iM 

[l/2 1/2)- 

Thus the two columns are equal, which means that it doesn't matter if ix jumps between j and 
j + 1, showing that v is continuous at xq. □ 
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8. Bounded operator fields 

In our investigation of Hessenberg reduction for matrix fields, our main problem has been 
that in the last few columns, as we approached the bottom-right of the matrix, we didn't have 
enough room anymore to continue the process, so we had to stop. This premature ending of the 
Hessenberg reduction process has been the main limiting factor to the strength of the results 
that we subsequently obtained, and one may wonder about ways to work around it. Of course, 
the premature stopping can't be completely avoided in general, because if we could obtain a 
complete Hessenberg form H^, that would in particular imply that we can completely separate 
eigenvalues, something that we already said is impossible in general over spaces of dimension 
3 and up. 

Another way of working around this problem is to add zeros to the right and to the bottom 
of our matrix field, so as to make more room. Going one step further one may add infinitely 
many rows and columns of zeros, making our matrix field a compact operator field; going one 
more step further, one may start right away with any bounded operator field. 

Throughout this section, let H = ^^(N*). Let B{H) denote the set of bounded operators on 
H. Let (ei)jgN* be its standard Hilbert basis. For x G B{H) and i,j £ N*, let 

be the (ij)-th matrix entry of x. 

Definition 8.1. A bounded operator x G B{H) is a Hessenberg operator if Xij = whenever 
i > j + 1 and Xij G M>o whenever i = j + 1. 

Definition 8.2. A Jacobi operator is a self-adjoint Hessenberg operator. 

Definition 8.3. Let H be a Hilbert space. A bounded operator x G B{H) is a cyclic operator 
if there exists a vector ^ £ H such that the linear span of the family {x^£,)k&^ is dense in H . 
The vector ^ is then called a cyclic vector of x. 

It is well-known (see [3j, Chapter 18) that the set of cyclic operators on H is not norm-dense 
in B{H), and that its complement is norm-dense. 

The following result is classical: 

Lemma 8.4. A bounded operator x G B{H) is cyclic if and only if it is of the form uhu* with 
u a unitary operator and h a Hessenberg operator. 

Proof. If h is Hessenberg then ei is a cyclic vector of it, hence uei is a cyclic vector of uhu*. 
Conversely, suppose that x is a cyclic operator. Let ^ be a cyclic vector of x. For k £ W, let 
fk = x^~^S,- Then {fk)keN* is a linearly independent family whose span is dense in H. Let 
{g^)keN* be the Hilbert basis of H obtained by applying the Gram-Schmidt orthonormalization 
process to (/fc)fceN*- Let u G B{H) be the unitary defined by letting u{ek) = g^ . It is then 
straightforward to check that u*xu is Hessenberg. □ 

Notice that it is a corollary of the previous lemma that any eigenspace of a cyclic operator 
has dimension 1. The converse is true for diagonalizable finite matrices. These facts suggest 
to think of cyclicity as a generalized notion of "multiplicity 1" . 

In this section, we prove the following theorem: 

Theorem 8.5. Let X be either a compact space or a finite- dimensional normal space. Let f 
be a strongly continuous map from X to B(H). For any e G C(X, M>o), there exist maps g and 
V from X to B{H) such that, letting h = v*gv, the following conditions hold: 
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• v{x) is an isometry for all x £ X . 

• h{x) is Hessenherg for all x £ X . 

• V, g, h are strongly continuous. 

• ll/(^) ~ 5(^)11 < ^(^) ^ G Moreover, f — g is norm- continuous, compact, and 
self- adjoint. 

We first need a few lemmas: 

Lemma 8.6. Let X be either a compact space or a finite- dimensional normal space. Let 
f G C{X,H) and e € C(X, M>o). It follows that there exists g £ C{X,H) such that for all 
xeX, \\f{x) - g{x)\\ < e{x) and g{x) ^ 0. 

Proof. The case when X is normal and finite-dimensional reduces to Lemma [1.6i Let us handle 
the case when X is compact. For /c G N* let be the set of all x G X such that 

(11) |(/(x),e,)| < e(x)/2 for alH ^ fe. 

The sequence {Uk)keN* is an open covering of X, and is increasing in the sense of inclusion. 
Since X is compact, it follows that there exists some k G N* such that Uk = X, so that 
inequality (llip holds for all x £ X. We may therefore obtain the wanted map g just by letting, 
for all X G X, 

g{x) = f{x) -^e{x)ek. □ 

The main technical step toward Theorem 18.51 is the following lemma: 

Lemma 8.7. Let X be either a compact space or a finite- dimensional normal space. Let f 
be a strongly continuous map from X to B{H). For any e G C(X, M>o), there exist sequences 
{9^)k€N* o.nd {u'^)k&i* of maps from X to B{H) such that, letting = u^g^u^* , the following 
properties are satisfied: 

• for all A; G N*, g^ and are strongly continuous maps from X to B{H) 

• for all /c G N*, is a norm- continuous map from X to U{H). 

• for all k G N*, g^^^ — g^ is norm- continuous, self-adjoint, compact, and \\g^~^^{x) — 
5^X2;) II < e{x)/2''. 

• for all k,l £ N* , if k ^ I then the k first rows of and agree. In other words, 
yi,j, k,l gW, ifi^k^l then u[j = u^j. 

• for all k,l N*, if 1 < k ^ I then the k — 1 first columns of and agree, and are 
in Hessenberg form. In other words: 

- Vi, j, k,l gW, if j <k^l then h\ - = h^-. 

- yi,j, k G N*, if j < k and i'^ j + 2 then h^ ■ = 0. 

- V«, j, k G N*, if j ^ k, j < n and i = j then h\ ■ / 0. 

Proof. We construct the sequences {g^)k&w-' and {u^)k(m* by induction. We set g^ = f, and we 
let be the constant map, u^{x) = id for all x G X. Now suppose that for some k G N*, the 
sequences u^, . . . , u'^ and g^ , . . . ,g'^ have been constructed, and let us then construct g^^^ and 
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Since the k — 1 first columns of are in Hessenberg form, that means that has the 
following form: 




In this diagram we represented the k — 1 first columns ending with zeros. Of course this only 
is an accurate representation if A; > 1. In the k = 1 case, one should simply imagine that the 
framed column is the first column. 



Let us do the same kind of "Householder reflection" argument that we already made in the 
proof of Theorem 13. 2[ In the above diagram, we framed the infinite block in the k-th column 
starting at row A; + 1. Call it b: 

oo 

Kx) = ^/ifc+i(x)ei. 

i=l 

Thus 6 is a map from X to H. Moreover, since is strongly continuous, b is continuous. Now 
let b' be the map obtained by taking off the first coefficient in b. In other words: 

oo 

(12) b'{x) = ^ht^,^,{x)e,. 

1=1 

By Lemma 18.61 ^ small perturbation ensures that b'{x) ^ for all x £ X. Moreover, this can 
be achieved by a self-adjoint perturbation on by applying the adjoint perturbation on the 
other side of the diagonal. So from now on we assume, without loss of generality, that h'' is 
such that the block b' inside it, as defined in equation ()12p . satisfies b'{x) ^ for all x £ X. 
This ensures that b{x) ^ and that 

c(^) •= Wr\u + ei / for all x G X 

The latter equation defines a continuous map c : X ^ H — {0}. Let vo{x) E U(//) denote the 
orthogonal refiection around the line generated by c(x). Notice that for all x € X, vo{x)b{x) 
belongs to the span of ei. Moreover, since c is continuous, vq is norm-continuous. Now define 
another unitary field v by letting, for all x £ X, 

It follows that u is still norm-continuous and that vh^v* satisfies all the wanted properties for 
/i'^"'"^. Let us define u'^'^^ := vvl^ , and g^~^^ = ■ That definition of g^^^ may look surprising 
but the truth is that we allowed ourselves above to modify "in place" , so of course it should 
be understood that g^ was then modified accordingly to preserve the equality = vl^g^u^* . 
At the next iteration, we have 

and we checked above that vh^v* has all the wanted properties for h^^^. □ 
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Lemma 8.8. Let X he a topological space. Let {f'^)k&i* be a sequence of strongly continuous 
maps from X to B[H). Suppose that there exists M G C(X, M>o) such that 

\\f^{x)\\ ^ M(x) for allx£ X and all k G N*. 

Suppose that for all /c E N* there exists Iq G N* such that for all I ^ Iq the first k columns of 
and of agree, in other words: 

flj{x) = fi]{x) for allxeX,l^j ^k,l^i. 

It follows that for all x G X, the sequence {f^{x))keN* converges strongly, and that, letting f{x) 
denote its strong limit, the map f thus defined on X is strongly continuous. 

Proof. For A; G N*, let : ^ denote the orthogonal projection onto the linear span of 
ei, . . . , Cfe. Let X G X. Let us check that {f^{x))ke^* converges strongly. Let ^ G if and 6 > 0. 
There exists G N* such that, letting ^' = p,,{^), we have - C|| < S/M{x). For ah / ^ k, 
we have f^x)^' = f^{x)^' , whence — /^(x)^'|| < 5 showing that the sequence (/^(x)^) 

is Cauchy. This shows that (/'(a^))/eN* converges strongly. The strong limit /(x) still satisfies 
||/(x)KM(x). 

Let us now prove that / is strongly continuous. Let xq G X. Let £ H. Let 5 > 0. There 
exists A: G N* such that, letting ^' = Pk{£,), we have 

By construction, for all I ^ k, we have f\x)^' = /'^(x)^' for all x G X. It follows that 
/(x)^' = /'^(x)^'. Since f^ is strongly continuous, this shows that the map a : x i— > /(x)^' 
is continuous on X. Let U be the set of all x G X such that ||a(x) — a(xo)|| < 5/2 and 
M(x) < 2M{xq). Thus U is an open neighborhood of xq of X. For all x £ U, we have 

\\f{x)C - /(xo)el ^ - fix)e\\ + \Hx) - a(xo)|| + ||/(xo)^' - /(xo)e||. 

But for all X G [/, 

ii/(x)c - /(x)c'K M(x)iie' - eiK V4. 

It follows that for all x G C/, 

\\f{x)C-fixom^S/A + 5/2 + 6/A^5. 

This proves that / is strongly continuous. □ 

Let us now finish the proof of our main theorem. 

Proof of Theorem \8. 5\ Apply Lemma [8.71 to obtain sequences {g^)k&i* and {u^)k&i* and let 
= u^g^u^*. By Lemma 18. 8| for all x G X, the sequences {u^*{x)) and {h^{x)) converge 
strongly, and, letting v{x) and h{x) denote their respective strong limits, the maps h thus 
defined on X are strongly continuous (we will prove below that h = v*gv). 

Notice that for all x G X, since t;(x) is an isometry because it is a strong limit of unitaries. 

Since 5'^(x)|| < e{x)/2^ and g^^^—g^ is self-adjoint, compact and norm-continuous 

for all k and cp = /, it follows that [g^ — f)k£N* converges locally uniformly to g — f and that 
g — f is self-adjoint, compact, norm-continuous, and ||5(x) — /(x)|| < e(x) for all x G X. This 
shows in particular that g is strongly continuous. 
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Let US show that h = v*gv. It is enough to show that (u g u ){x) converges weakly to 
{v*gv){x) for all x & X. Let ^,r] e H. We have 

((tx^/u'=*)(x)e,r?) = ((5^'=*)(x)^,(tx'=*)(x)r?). 

As the sequences (^f^) and (u^*) converge strongly to g and v respectively, and the strong 
topology makes multiplication jointly continuous on bounded sets, the sequence {g^u^*) also 
converges strongly to gv, and we obtain 

{iu''g''u''*)ix)^,r]) ^ {igv)ix)^,vix)r]) = {{v* gv)ix)^,r]) . 

This proves that h = v*gv. □ 
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